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ABSTRACT 

We  report  on  an  outburst  of  the  high  mass  X-ray  binary  4U  0115-^634  with  a  pulse  period  of  3.6  s  in  2008  March/April  as  observed 
with  RXTE  and  INTEGRAL.  During  the  outburst  the  neutron  star’s  luminosity  varied  by  a  factor  of  10  in  the  3-50  keV  band.  In 
agreement  with  earlier  work  we  find  evidence  of  five  cyclotron  resonance  scattering  features  at  ~10.7,  21.8,  35.5,  46.7,  and  59.7  keV. 
Previous  work  had  found  an  anticorrelation  between  the  fundamental  cyclotron  line  energy  and  the  X-ray  fiux.  We  show  that  this 
apparent  anticorrelation  is  probably  due  to  the  unphysical  interplay  of  parameters  of  the  cyclotron  line  with  the  continuum  models 
used  previously,  e.g.,  the  negative  and  positive  exponent  power  law  (NPEX).  For  this  model,  we  show  that  cyclotron  line  modeling 
erroneously  leads  to  describing  part  of  the  exponential  cutoff  and  the  continuum  variability,  and  not  the  cyclotron  lines.  When  the 
X-ray  continuum  is  modeled  with  a  simple  exponentially  cutoff  power  law  modified  by  a  Gaussian  emission  feature  around  10  keV, 
the  correlation  between  the  line  energy  and  the  fiux  vanishes,  and  the  line  parameters  remain  virtually  constant  over  the  outburst. 

We  therefore  conclude  that  the  previously  reported  anticorrelation  is  an  artifact  of  the  assumptions  adopted  in  the  modeling  of  the 
continuum. 

Key  words.  X-rays:  binaries  -  pulsars:  individual  4U  0115-^634  -  magnetic  fields 


1.  Introduction 

In  all  accreting  X-ray  pulsars  the  accretion  geometry  is  strongly 
affected  by  the  magnetic  field  of  the  neutron  star,  which  has  val¬ 
ues  of  a  few  10^^  G  at  the  magnetic  poles.  The  accreted  matter 
is  forced  to  follow  the  magnetic  field  lines  from  the  Alfven  ra¬ 
dius  inward  to  the  accretion  columns  located  at  the  magnetic 
poles  of  the  neutron  star.  The  shape  of  these  columns  remains 
an  open  question.  The  simplest  models  assume  cylinder  or  slab¬ 
like  geometries,  although  more  complicated  shapes  have  been 
discussed  (Meszaros  1984;  Becker  et  al.  2012).  On  each  mag¬ 
netic  pole  of  the  neutron  star,  a  hot  spot  is  created,  emitting  ther¬ 
mal  photons  into  the  accretion  column  (Davidson  1973).  These 
photons  are  upscattered  in  energy  by  inverse  Compton  scatter¬ 
ing  within  the  accretion  column,  forming  the  hard  X-ray  spec¬ 
trum  (Becker  &  Wolff  2007;  Ferrigno  et  al.  2009).  Depending 
on  the  source  luminosity,  and  thus  mass  accretion  rate,  M,  and 
deceleration  mechanism,  a  radiation-dominated  shock  can  form 


above  the  surface  of  the  neutron  star.  See  Becker  et  al.  (2012)  for 
a  discussion  of  the  different  modes  of  accretion  in  the  accretion 
column. 

Owing  to  the  strong  magnetic  field  of  accreting  X-ray  pul¬ 
sars,  cyclotron  resonance  scattering  features  (cyclotron  lines; 
CRSFs)  are  observable  in  many  X-ray  pulsars.  These  features 
have  been  reported  for  about  20  X-ray  pulsars  (Caballero  & 
Wilms  2011).  CRSFs  are  produced  by  photons  generated  in  the 
accretion  column  of  the  neutron  star  interacting  with  electrons 
(Schonherr  et  al.  2007;  Meszaros  1992,  and  references  therein). 
The  motion  of  these  electrons  perpendicular  to  the  5-field  is 
quantized  into  Landau  levels  with  energy  differences  of  about 


Rhotons  inelastically  scattering  off  these  electrons  will  imprint 
absorption-line  like  features  onto  the  spectral  continuum. 
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Measuring  the  energy  of  the  CRSF  yields  information  about 
the  field  strength  in  the  line  forming  region.  Since  the  5-field 
of  a  neutron  star  can  be  assumed  to  be  dipole-like  to  first  order, 
in  many  models  the  centroid  energy  of  CRSFs  depends  on  the 
characteristic  emission  height,  i.e.,  the  altitude  of  the  line  form¬ 
ing  region  (Basko  &  Sunyaev  1976;  Burnard  et  al.  1991;  Mihara 
et  al.  2004).  Since  this  altitude  depends  on  the  complex  rela¬ 
tion  between  the  deceleration  of  the  accreted  material  and  the 
mass  accretion  rate,  it  is  expected  to  depend  on  the  source  lu¬ 
minosity.  Such  a  dependence  has  indeed  been  observed.  Three 
types  of  behavior  have  been  observed.  One  type  shows  a  pos¬ 
itive  correlation  between  the  X-ray  luminosity  Lx  and  the  cy¬ 
clotron  line  energy  E^yc  (e.g.,  Her  X-1,  Staubert  et  al.  2007). 
The  other  type  shows  a  negative  correlation  between  Lx  and  E^yc 
(e.g.,  V0332-h53,  Nakajima  et  al.  2010,  and  references  therein). 
Finally,  there  are  also  cyclotron  sources  known  with  a  constant 
cyclotron  line  energy  (e.g.,  A0535-r26,  Caballero  et  al.  2007). 

As  suggested  by  Staubert  et  al.  (2007)  and  Klochkov  et  al. 
(2011)  and  then  quantified  by  Becker  et  al.  (2012),  these  differ¬ 
ent  relations  between  Lx  and  E^yc  could  be  due  to  the  different 
deceleration  mechanisms  for  the  accreted  material.  Above  a  cer¬ 
tain  critical  luminosity.  Lent,  the  accreted  material  is  decelerated 
mainly  by  the  radiation  pressure.  In  this  situation,  an  increase  in 
luminosity  stops  the  material  earlier,  leading  to  the  line  forming 
region  being  in  a  region  of  lower  5- field  strength  and  yielding 
a  negative  correlation  between  E^yc  and  Lx-  For  very  low  lumi¬ 
nosities,  the  deceleration  is  dominated  by  the  gas  pressure.  In 
this  regime  higher  mass  accretion  rates  will  compress  the  accre¬ 
tion  mound  and  thus  lead  to  a  positive  correlation  between  the 
luminosity  and  the  cyclotron  line  energy.  Finally,  at  intermediate 
luminosities  one  expects  a  regime  where  the  accretion  column 
is  fairly  stable  and  the  deceleration  is  dominated  by  Coulomb 
interactions. 

4U  01 15-r634  (Giacconi  et  al.  1972)  is  one  of  the  X-ray  pul¬ 
sars  for  which  CRSFs  have  been  studied  in  great  detail  (see, 
e.g.,  Wheaton  et  al.  1979;  White  et  al.  1983;  Nagase  et  al.  1991; 
Heindl  et  al.  1999;  Santangelo  et  al.  1999;  Mihara  et  al.  2004; 
Nakajima  et  al.  2006;  Tsygankov  et  al.  2007;  Ferrigno  et  al. 
2009,  2011).  The  source  shows  3.61  s  pulsation  and  has  an  or¬ 
bital  period  of  -24.3  d  (Cominsky  et  al.  1978;  Rappaport  et  al. 
1978).  The  optical  counterpart  is  the  09e  star  V635  Cas  (Johns 
et  al.  1978;  Unger  et  al.  1998). 

In  previous  outbursts,  CRSFs  have  been  detected  up  to  the 
fifth  harmonic  (Heindl  et  al.  2000).  This  is  the  largest  number 
of  detected  CRSFs  in  an  accreting  X-ray  pulsar.  In  later  work, 
a  negative  correlation  between  Lx  and  E^yc  was  found  (e.g., 
Nakajima  et  al.  2006;  Tsygankov  et  al.  2007 ;  Muller  et  al.  2010), 
meaning  that  the  source  should  be  located  in  the  supercritical  lu¬ 
minosity  regime.  Other  authors  even  found  a  discontinuous  be¬ 
havior  of  the  fundamental  line  energy  and  explained  this  by  a 
sudden  change  of  the  accretion  conditions  (Li  et  al.  2012). 

In  this  work  we  revisit  the  behavior  of  the  CRSFs  in 
4U  01 15-r634.  Focusing  on  data  taken  during  a  giant  outburst  in 
2008  March/ April  (Fig.  1).  We  show  that  the  previously  claimed 
relation  between  the  source  flux  Lx  and  E^yc  strongly  depends 
on  the  choice  of  the  underlying  continuum  model.  In  Sect.  2  we 
summarize  the  observations.  In  Sect.  3  we  describe  the  differ¬ 
ent  continuum  models  and  the  flux  and  time  dependence  of  the 
continuum  and  cyclotron  line  parameters.  We  argue  that  the  cy¬ 
clotron  line  is  only  modeled  correctly  when  using  an  exponen¬ 
tially  cutoff  power  law  continuum  with  a  strong  emission  fea¬ 
ture  around  10  keV.  Furthermore,  alternative  continua  models, 
which  have  been  used  in  some  earlier  investigations,  result  in  in¬ 
correct  cyclotron  line  parameters  which  can  cause  the  apparent 


relation  between  source  flux  Lx  and  fundamental  cyclotron  line 
energy  Lq.  We  summarize  this  paper  in  Sect.  5. 

2.  Observations  and  data  reduction 

The  outburst  of  4U  0115-r634  started  2008  March  12,  after  a 
phase  of  almost  four  years  of  quiescence,  when  Swiftj^NI  de¬ 
tected  a  significant  increase  of  the  X-ray  flux  (Krimm  et  al. 
2008).  The  outburst  lasted  about  40  days  and  exceeded  a  flux 
of  -200mCrab  in  the  2-12  keV  RXTE/ASM  band  (Fig.  1).  The 
whole  outburst  was  monitored  with  pointed  observations  by  the 
Rossi  X-ray  Timing  Explorer  (RXTE,  Bradt  et  al.  1993)  and 
several  pointed  observations  by  the  INTErnational  Gamma-Ray 
Astrophysics  Laboratory  (INTEGRAL,  Winkler  et  al.  2003). 
Here,  we  use  data  from  RXTE’s  Proportional  Counter  Array 
(PGA,  Jahoda  et  al.  2006)  and  the  High  Energy  X-Ray  Timing 
Experiment  (HEXTE,  Rothschild  et  al.  1998),  as  well  as  the 
INTEGRAL  Soft  Gamma-Ray  Imager  (ISGRI,  Lebrun  et  al. 
2003).  Data  were  reduced  with  our  standard  analysis  pipelines, 
based  on  HEASOET  (v.  6.11)  and  INTEGRAL  OSA  v.  9.0. 
Data  modeling  was  performed  with  the  newest  version  of 
the  Interactive  Spectral  Interpretation  System  (ISIS,  Houck  & 
Denicola  2000). 

The  PCA  consisted  of  five  proportional  counter  units 
(PCUs),  which  were  sensitive  between  2  and  60  keV.  Since 
PCU  2  is  known  to  have  been  the  best  calibrated  one  (Jahoda 
et  al.  2006),  we  only  extracted  spectra  from  the  top  layer  of  this 
PCU  in  the  standard2f  mode.  The  hard  X-ray  spectrum  from 
15-250  keV  was  monitored  with  HEXTE.  HEXTE  consisted  of 
two  independent  clusters,  A  and  B,  both  arrays  of  four  phoswich 
scintillation  detectors.  The  “rocking  mechanism”  of  these  clus¬ 
ters  allowed  a  near  realtime  measurement  of  the  background. 
Since  this  mechanism  failed  for  cluster  A  before  our  observa¬ 
tions,  we  only  used  data  from  cluster  B.  The  other  instrument 
for  the  hard  X-rays  was  the  pixelized  CdTe  detector  ISGRI  be¬ 
hind  a  Coded  Mask  on  board  the  INTEGRAL  satellite.  It  covers 
the  energy  range  between  -18  keV  and  IMeV.  We  extracted 
all  available  ISGRI  spectra,  for  which  4U  0115-1-634  was  less 
than  10°  off-axis.  Table  1  contains  an  observation  log  (see  also 

Fig.  1). 

Eor  the  spectral  analysis  we  used  data  between  3  and  50  keV 
for  PCA,  and  between  20  and  100  keV  for  HEXTE  and  for 
ISGRI.  A  systematic  error  has  been  added  in  quadrature  to 
the  PCA  and  ISGRI  data,  using  the  canonical  values  of  0.5% 
and  2.0%  (Jahoda  et  al.  2006,  and  IBIS  Analysis  User  Manual^), 
respectively.  To  improve  the  signal  to  noise  ratio  of  the  spectra, 
we  averaged  the  data  over  14  data  blocks  in  time.  These  data 
blocks  are  defined  in  the  fifth  column  of  Table  1  and  shown  in 
Pig.  1.  We  omitted  one  ISGRI  dataset  from  our  analysis  (rev¬ 
olution  number  677),  because  no  simultaneous  RXTE  data  are 
available  at  that  time  and  the  source  was  almost  in  the  off  state 
during  that  observation. 

3.  Spectral  analysis 

3.1.  Pulsar  continuum  models 

A  large  problem  when  analyzing  X-ray  spectra  from  X-ray  bi¬ 
nary  pulsars  is  to  find  a  good  description  of  the  broadband  con¬ 
tinuum.  In  most  cases,  semi-physical  models,  which  are  some 
variant  of  a  power  law  with  a  high  energy  exponential  cutoff 

^  http : //WWW . isdc . unige . ch/integral/analysis# 
Documentation 
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Table  1.  Log  of  observations  for  the  2008  outburst  of  4U  0115-^634. 


ID« 

Start  date 

MJD 

^exp^  [s] 

e^ 

01-01-00 

2008  Mar.  20 

54545.27-54  545.43 

8336/2413 

1 

01-01-01 

2008  Mar.  20 

54545.47-54  545.49 

1648/382 

1 

01-02-04 

2008  Mar.  21 

54546.97-54  547.20 

12448/3938 

2 

01-02-00 

2008  Mar.  23 

54548.07-54  548.25 

9904/3082 

2 

01-02-01 

2008  Mar.  24 

54549.19-54  549.30 

6336/1948 

3 

01-02-02 

2008  Mar.  24 

54549.97-54  550.14 

8832/3136 

3 

01-02-03 

2008  Mar.  25 

54550.24-54  550.28 

2384/934 

3 

01-02-05 

2008  Mar.  26 

54551.08-54  551.17 

4096/1679 

4 

01-02-06 

2008  Mar.  26 

54551.28-54  551.31 

752/625 

4 

01-03-02 

2008  Mar.  29 

54554.23-54  554.38 

576/2112 

5 

01-03-01 

2008  Mar.  31 

54556.14-54  556.16 

768/473 

6 

01-03-00 

2008  Mar.  31 

54556.25-54  556.37 

768/1558 

6 

01-03-03 

2008  Apr.  02 

54558.93-54  559.10 

5456/3104 

7 

01-04-00 

2008  Apr.  04 

54560.04-54  560.22 

4240/2998 

7 

01-04-01 

2008  Apr.  06 

54562.01-54  562.18 

5888/2908 

8 

01-04-02 

2008  Apr.  06 

54562.95-54  563.09 

5552/2317 

8 

01-04-03 

2008  Apr.  07 

54563.99-54  564.07 

3856/1548 

9 

01-04-04 

2008  Apr.  08 

54564.82-54  564.98 

7808/2853 

9 

01-04-05 

2008  Apr.  09 

54565.86-54  566.10 

10  896/4101 

10 

01-04-06 

2008  Apr.  10 

54566.98-54  567.00 

1648/699 

10 

01-05-00 

2008  Apr.  11 

54567.00-54  567.01 

1056/381 

11 

01-05-01 

2008  Apr.  11 

54567.96-54  568.13 

9632/3176 

11 

16-01-00 

2008  Apr.  13 

54569.98-54  570.16 

10  192/3043 

12 

16-01-01 

2008  Apr.  15 

54571.94-54  572.12 

10  368/3162 

12 

16-02-00 

2008  Apr.  18 

54574.04-54  574.15 

7184/1924 

13 

16-02-01 

2008  Apr.  22 

54578.36-54  578.51 

6480/1679 

14 

664 

2008  Mar.  21 

54546.48-54  547.68 

55  981 

2 

667 

2008  Mar.  30 

54555.46-54  556.04 

28  279 

6 

668 

2008  Apr.  02 

54558.45-54  559.05 

27  655 

7 

669 

2008  Apr.  07 

54563.38-54  564.02 

28  495 

9 

670 

2008  Apr.  10 

54566.41-54  567.01 

30990 

10 

673 

2008  Apr.  17 

54573.90-54  574.50 

31360 

13 

675 

2008  Apr.  24 

54580.62-54  581.97 

70902 

14 

677 

2008  Apr.  29 

54585.69-54  586.91 

63  662 

- 

Notes.  The  first  part  of  the  table  contains  the  RXTE  observations,  the 
second  part  the  ISGRI  data.  For  RXTE,  the  first  column  contains 
the  number  of  the  Obs-ID  after  “93032-”.  For  ISGRI,  it  contains  the 
revolution  number.  For  RXTE,  the  two  numbers  correspond  to  the 
PCA  and  HEXTE  exposures,  respectively.  Epoch  for  data  grouping, 
see  text  for  details. 


(Kreykenbohm  et  al.  2002,  and  references  therein),  are  used  to 
describe  the  spectra.  These  models  can  be  justified  by  consider¬ 
ing  the  Comptonization  of  soft  photons  in  an  accretion  column 
(Becker  &  Wolff  2007,  and  references  therein).  The  most  basic 
exponentially  cutoff  power  law  model,  called  Cutoff  PL  in  ISIS 
and  XSPEC  (Arnaud  1996),  has  the  form 

CutoffPL(£’)  cx  E~^  exp  (-E’/E’foid) ,  (2) 


1998  2000  2002  2004  2006  2008  2010 


1997  1999  2001  2003  2005  2007  2009 


Fig.  1.  RXTE/ASM  lightcurve  of  4U  0115-^634  (2-12  keV).  The  inset 
shows  a  close-up  view  of  this  lightcurve,  showing  the  2008  outburst 
in  detail.  Arrows  indicate  the  times  of  pointed  RXTE  and  INTEGRAL 
observations.  Numbers  indicate  the  observation  epochs  used  in  the  data 
analysis  (see  also  Table  1). 


where  for  typical  X-ray  pulsars  E'break  ~  10  keV  and  which  has 
an  unphysical  “kink”  at  E’break,  which  can  result  in  line-like  resid¬ 
uals  around  this  energy  (Kretschmar  et  al.  1997;  Kreykenbohm 
et  al.  1999).  For  this  reason  Klochkov  et  al.  (2008)  and  Ferrigno 
et  al.  (2009)  smooth  the  continuum  around  this  break,  using  a 
model  of  the  form 


PLINT(£)  oc 


E-^ 

AE^  +  BE^  +  CE  +  D 
exp  (-E’/E’foid) 


for  E  <  Ebreak  “  AE 
for  |E  -  EbreakI  <  AE 
for  E  >  Ebreak  +  AE 


(4) 


where  the  coefficients  A,  E,  C,  and  D  are  chosen  such  that  the 
continuum  and  its  first  derivative  are  continuous  at  E  =  Ebreak  ± 
AE  and  where  AE  0  but  sufficiently  small. 

Other  common  choices  to  avoid  the  break  are  the  so  called 
negative  and  positive  power-law  exponential  (NPEX,  see,  e.g., 
Makishima  et  al.  1999,  and  references  therein),  a  continuum 
model  with  five  free  parameters  of  the  form 

NPEX(£:)  (X  (£"^‘  +  exp  {-EIEma)  ,  (5) 

with  photon  indices  >  0^,  and  the  “Fermi  Dirac  cutoff” 

(FDCUT,  see  Tanaka  1986) 


with  the  photon  index,  F,  and  the  folding  energy,  Efoid  (see,  e.g., 
Schonherr  et  al.  2007).  However,  this  simple  model  is  often, 
but  not  always,  insufficient  to  describe  the  continua  of  observed 
X-ray  pulsars,  as  the  curvature  in  the  continuum  is  often  seen  to 
start  at  different  energies.  To  obtain  a  good  spectral  fit,  therefore, 
models  with  more  free  fit  parameters  and  more  complex  shapes 
are  sometimes  required.  The  most  simple  of  these  models  is  the 
power  law  with  high  energy  cutoff 


PLCUT(E)  oc 


E“^ 

E“^  exp  ((Ebreak  “  F’)/^fold) 


for  E  <  Ebreak 
for  E  >  Ebreak 


(3) 


FDCUT(E)  oc  E“^ 


1 

1  -k  exp  ((E  -  Ebreak)/F’fold) 


(6) 


Common  to  all  of  these  empirical  continua  is  that  despite  the 
smoother  transition  around  Ebreak  they  often  need  to  be  modified 
further  by  a  broad  excess  or  depression  around  10  keV.  This  so- 
called  “10  keV  feature”  has  been  observed  in  many  X-ray  pulsar 
spectra,  either  in  absorption  or  in  emission,  although  its  origin 
is  still  not  clear  (Coburn  2001).  This  feature  is  usually  modeled 

^  Often  one  fixes  T2  -  2,  see,  e.g.,  Nakajima  et  al.  (2006). 
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as  an  additive  Gaussian  line,  with  the  centroid  energy  Eq,  the 
width  (Tg,  and  the  flux  Aq. 

In  addition  to  those  empirical  models,  based  on  the  the¬ 
oretical  description  of  bulk  and  thermal  Comptonization  in 
magnetized  accretion  columns  of  Becker  &  Wolff  (2007,  and  ref¬ 
erences  therein),  Ferrigno  et  al.  (2009)  have  developed  a  contin¬ 
uum  model  depending  on  physical  parameters  like  the  column 
radius,  the  accretion  rate,  the  electron  temperature,  the  5-fleld 
strength.  A  hybrid  model  which  combines  the  broadband  contin¬ 
uum  with  CRSFs  is  currently  in  preparation  (Schwarm  et  al.,  in 
prep.).  These  physics-based  continuum  models  are  still  in  active 
development  and,  because  of  the  large  number  of  free  parame¬ 
ters,  only  suited  for  very  high  signal  to  noise  data.  For  this  reason 
most  observational  work  uses  one  of  the  continua  described  by 
Eqs.  (2)-(6). 

For  sources  with  cyclotron  lines,  the  continuum  model  is 
modified  with  a  description  of  the  cyclotron  line.  Since  phys¬ 
ically  self  consistent  line  models  such  as  those  described  by 
Schonherr  et  al.  (2007)  are  still  not  available  for  all  sources, 
empirical  line  models  are  generally  used.  The  two  most  com¬ 
mon  CRSF  models  are  those  using  a  Gaussian  or  a  pseudo- 
Lorentzian  profile  (see  discussion  by  Nakajima  et  al.  2010). 
Here,  we  model  the  CRSFs  with  a  pseudo-Lorentzian  optical 
depth  profile  (Mihara  et  al.  1990),  i.e.,  we  multiply  the  contin¬ 
uum  model  with 


CYCLABS(^)  =  exp 


T{WEIE,y,f  \ 

(£-£cyc)2  +  w7’ 


(7) 


with  the  centroid  energy,  Ecyc,  the  width  of  the  feature,  W,  and 
the  optical  depth  of  the  line,  r.  This  approach  allows  the  com¬ 
parison  of  the  line  parameters  with  results  from  previous  papers, 
in  which  also  the  CYCLABS  model  has  been  used  (e.g.,  Nakajima 
et  al.  2006;  Tsygankov  et  al.  2007 ;  Li  et  al.  2012).  In  the  remain¬ 
der  of  the  paper,  we  label  the  cyclotron  line  parameters  with  the 
number  of  respective  harmonic,  where  0  corresponds  to  the  fun¬ 
damental  line.  As  common  for  CRSFs,  W  and  r  are  strongly  cor¬ 
related  (Coburn  2001).  We  therefore  set  a  lower  limit  of  0.5  keV 
for  the  width,  which  is  comparable  to  the  typical  PCA  resolution. 

In  order  to  describe  the  data  of  4U  01 15-h634,  in  addition  to 
several  of  the  continuum  models  discussed  above,  we  take  in¬ 
terstellar  absorption  into  account,  modeling  it  with  an  updated 
version  of  TBabs^,  using  abundances  by  Wilms  et  al.  (2000)  and 
cross  sections  by  Vemer  &  Yakovlev  (1995).  We  furthermore  ac¬ 
count  for  an  intrinsic  Fe  Ka  fluorescence  line,  which  we  describe 
with  an  additive  narrow  Gaussian  emission  line  with  frozen  cen¬ 
troid  energy  at  6.4  keV,  a  width  of  lO""^  keV,  and  a  flux  Ape  and 
equivalent  width  Wpe-  Assuming  a  narrow  line  is  justified  be¬ 
cause  in  systems  like  4U  0115-r63  Fe  Ka  is  expected  to  have  a 
width  on  the  order  of  <0.5  keV  (see,  e.g.,  Torrejon  et  al.  2010), 
which  cannot  be  resolved  by  PCA. 

In  order  to  take  into  account  that  the  flux  normalization  of  the 
different  instruments  used  is  not  perfectly  known  we  introduce 
cross  calibration  constants  chexte  and  cisgri  as  free  fit  param¬ 
eters.  These  constants  also  account  for  flux  differences  of  the 
source  between  the  observations  which  were  not  fully  simulta¬ 
neous.  The  PCA  background  was  estimated  using  the  standard 
models  provided  by  GSFC.  These  estimates  often  show  a  flux 
deviation  from  the  proper  background  by  a  few  percent.  This  de¬ 
viation  was  accounted  for  by  multiplying  the  background  flux 
with  another  fit  parameter,  cp  (see  also  Rothschild  et  al.  201 1). 

^  See  http://pulsar.sternwarte.uni-erlangen.de/wilms/ 
research/tbabs/ 


energy  [keV] 

Fig.  2.  Deviations  between  the  data  and  the  best  fit  with  Cutof fPL  and 
10  keV  feature  for  the  PCA  (circles),  HEXTE  (crosses),  and  ISGRI 
(filled  diamonds)  in  units  of  standard  deviations  cr.  Eor  these  fits,  no 
cyclotron  lines  have  been  taken  into  account.  The  numbers  on  the  right 
side  indicate  the  respective  epoch. 


3.2.  Fitting  strategy 

Since  the  Cut  off  PL  is  the  easiest  continuum  model  for  X-ray 
pulsars,  we  used  this  model  for  our  fits.  Furthermore,  the  10  keV 
feature  in  emission  is  required  to  get  a  good  description  of  the 
data.  For  each  epoch,  the  fits  have  been  simultaneously  per¬ 
formed  for  all  available  instruments.  Figure  2  shows  the  residu¬ 
als  of  the  best  fit  using  this  continuum  without  taking  any  CRSFs 
into  account.  These  residuals  clearly  indicate  the  need  of  the 
cyclotron  lines  in  this  model  at  ~11,  22,  and  sometimes  above 
30  keV.  As  the  quality  of  individual  data  sets  is  strongly  variable 
over  the  outburst,  in  an  initial  fit  run  a  variable  number  of  CRSFs 
were  added  to  individual  fits  until  the  reduced  ~  1-  To 

avoid  continuum  modeling  with  the  absorption  features,  we  fixed 
the  widths  of  the  CRSFs  of  the  second  and  higher  harmonics  to 
typical  values,  i.e.,  4  keV  (Ferrigno  et  al.  2009). 

In  order  to  check  whether  the  number  of  different  cyclotron 
features  in  the  model  affects  the  fit  results  of  the  continuum 
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Nil  [10^2  cm  2]  Efoid  [keV] 


0  2  4  10  12 


0  2  4  10  11  12  13 


Nil  [10^2  cm-2]  Eim  [keV] 

Fig.  3.  a)  c)  and  d)  Confldence  contours  for  the  folding  energy,  the  pho¬ 
ton  index  and  the  hydrogen  column  density  for  epoch  6.  b)  Confldence 
contours  between  the  Ffoid  and  F  with  A^h  frozen  to  1.7  X  10^^  cm“^. 
Contour  lines  correspond  to  the  68.3%,  90%,  and  99%  level.  Color  indi¬ 
cates  Ax^  with  respect  to  the  best  fit  value,  with  the  color  scale  running 
from  orange  (low  Ax^)  to  dark  blue  (large  Ax^). 


parameters  and  the  fundamental  and  first  harmonic  CRSF,  in  a 
second  run  all  spectral  fits  included  only  the  fundamental  line 
and  its  first  and  second  harmonics.  For  epochs,  where  only  the 
first  harmonic  was  necessary  to  obtain  a  good  fit  in  our  initial 
run,  we  added  the  second  harmonic  with  fixed  parameters  dur¬ 
ing  the  second  run,  holding  its  parameters  at  the  mean  values  of 
the  results  from  the  other  epochs. 

Since  initial  fits  with  the  Cutoff  PL  continuum  showed  large 
uncertainties  for  the  hydrogen  column  density,  A^h,  without  a  sig¬ 
nificant  variability  between  different  epochs,  we  held  A^h  fixed  at 
the  mean  value  A^h  =  1-7  x  10^^  cm“^  to  avoid  unphysical  corre¬ 
lations  between  these  fit  parameters  (Fig.  3). 

The  two  fit  runs  show  that  almost  all  fit  parameters  are  equal 
within  their  respective  uncertainties  for  both  runs.  Only  £’foid  and 
the  line  parameters  of  the  second  harmonic  are  slightly  different 
in  epoch  6  and  7.  These  differences  can  be  explained  by  the  fact 
that  for  these  two  cases  the  second  run  has  to  account  for  the 
equivalent  of  the  four  harmonics  required  in  the  first  run.  The 
differences  affect  mainly  the  higher  energies  and  are  compen¬ 
sated  by  a  change  in  the  folding  energy  and  the  closest  cyclotron 
line,  i.e.,  the  second  harmonic.  The  remainder  of  this  paper  is 
based  on  the  results  from  the  first  fit  run. 

An  example  for  the  best  fit  using  the  Cut  off  PL  continuum 
with  10  keV  feature  and  five  cyclotron  features  (first  fit  run, 
epoch  6)  is  shown  together  with  all  model  components  in  Fig.  4. 
This  spectrum  was  recorded  by  all  three  instruments  during  the 
maximum  of  the  outburst,  i.e.,  epoch  6.  All  fit  results  from  runs 
using  the  Cutoff  PL  continuum  are  summarized  in  Table  2. 


4.  Spectral  evolution:  No  anticorrelation 
of  the  cyciotron  iine  energy  with  fiux 

Figures  5-7  show  the  evolution  of  the  continuum  and  cyclotron 
line  parameters  as  well  as  the  parameters  for  the  Fe  Ka  line  for 


energy  [keV] 

Fig.  4.  a)  Folded  spectrum  from  epoch  6  together  with  best  fit  (red) 
and  all  model  components,  i.e.,  CutoffPL,  10  keV  feature,  Fe  Ka  flu¬ 
orescence  line,  and  five  cyclotron  features,  b)  Residuals  of  the  best  fit, 
shown  as  the  difference  between  the  data  and  the  model  normalized  to 
the  uncertainty  of  each  data  point.  Dark  blue  data  points  correspond  to 
PCA,  cyan  points  to  HEXTE,  and  orange  points  to  ISGRI. 


the  first  fit  run  with  the  CutoffPL  continuum  and  10  keV  fea¬ 
ture,  and  fixed  Ah  (see  also  Table  2).  As  already  noted  in  previ¬ 
ous  analyses  for  4U  0115-r634,  the  continuum  varies  over  the 
outburst  (e.g.,  Li  et  al.  2012).  During  the  brightest  phases  of 
the  outburst  the  spectrum  is  the  hardest  one.  As  shown  in  Fig.  3, 
the  variation  of  F  is  significantly  larger  than  the  very  slight  cor¬ 
relation  seen  between  F  and  E’foid-  As  expected,  the  source’s  in¬ 
trinsic  iron  Ka  line  is  significantly  present  during  all  epochs. 
As  also  observed  in  other  X-ray  transient  pulsars  (Inoue  1985) 
and  expected  for  a  line  originating  from  fluorescence,  the  line 
fiux  is  positively  correlated  to  the  X-ray  fiux  from  the  source. 
The  equivalent  width  stays  constant  within  the  uncertainties  over 
the  outburst.  The  10  keV  feature  varies  in  both,  the  centroid  en¬ 
ergy  Eq  and  its  fiux  Aq.  Its  energy  seems  to  be  slightly  anticor¬ 
related  with  the  source  flux,  while  the  relative  flux  in  the  feature, 
i.e.,  the  ratio  between  the  flux  in  the  10  keV  feature  and  the  total 
flux,  remains  constant.  Its  width,  ctq,  remains  relatively  constant 
between  3.0  keV  and  3.5  keV. 

Very  different  from  virtually  all  previous  analyses,  however, 
the  cyclotron  line  is  seen  to  be  only  slightly  varying  over  the  out¬ 
burst.  The  colored  data  points  in  Figs.  6  and  7  show  the  centroid 
energy,  Eq,  the  width,  Wq,  and  the  optical  depth,  tq,  of  the  fun¬ 
damental  cyclotron  line  against  time  and  flux  for  the  CutoffPL 
continuum,  modified  with  the  10  keV  feature.  In  these  fits,  the 
width  Wo  remains  constant  throughout  the  outburst,  and  there  is 
only  a  slight  (~3(t)  indication  that  the  optical  depth  tq  was  some¬ 
what  shallower  during  the  initial  phases  of  the  outburst,  before 
~MJD  54552,  although  it  cannot  be  excluded  that  some  of  this 
effect  is  due  to  a  correlation  with  the  ‘TO  keV  feature”,  which 
during  this  time  of  the  outburst  peaks  close  to  the  cyclotron 
line  at  around  8  keV.  Most  importantly,  however,  as  shown  in 
Fig.  6  the  energy  of  the  fundamental  cyclotron  line  is  only  very 
slightly  variable:  during  early  phases  of  the  outburst  the  line 
is  around  11  keV,  moving  towards  10  keV  as  the  outburst  pro¬ 
gressed.  This  result  is  contrary  to  all  earlier  studies  of  outbursts 
of  4U  0115-r634,  where  strong  changes  in  the  line  energy  were 
seen  (e.g.,  Li  et  al.  2012;  Tsygankov  et  al.  2007 ;  Nakajima  et  al. 
2006;  Mihara  et  al.  2004). 

Where  does  this  difference  in  line  behavior  come  from? 
As  discussed  in  Sect.  3.1  a  large  variety  of  continuum  models 
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Table  2.  Continuum  parameters  of  the  time  resolved  spectral  analysis. 


ep^ 

T  b 

F 

f^fold 

Ape^ 

Wpe 

Eg 

Ag^ 

O-Q 

% 

ChEXTE 

ClSGRl 

+L/d.o-f- 

[keV] 

[10-3] 

[eV] 

[keV] 

[keV] 

1 

0  59+0-04 

10  7+0-5 

-0.5 

9  c+0.6 
^•'^-0.6 

55+12 

JJ-12 

Q  9+0.3 

O.J-02 

0-18!?:?? 

9  4+0.1 
-0.3 

1  00+^-^^ 
4-44^-0.02 

0  73+0-02 
'  ^-0.02 

- 

1.13/50 

2 

c  a  c +0.02 
^'•^^-0.02 

-0.03 

1 1  5+0-4 

^^•^-0.4 

9  q+0.7 
^•^-0.7 

53+^^ 

8-01? 

0.19!?:?? 

9  9+0.1 
•+-•+-0.1 

0  96+^-^^ 

44-^'+_qqj 

0  74+0-01 

44-  '^-0.01 

0  0 
d  d 
+  1 

00 

0 

1.12/77 

3 

7.48t“ 

0  54+0-03 
-0.03 

11.61:? 

9  7+0.9 
^’'-0.9 

46+10 
^  -10 

7  9+0.2 
' -0-0.2 

0  34+0-02 
-0.01 

3  6+^-^ 
•+•'+-0.1 

0  94+0-02 
44- ^^-0.02 

0.75!?:?! 

- 

1.32/48 

4 

8.66:«:“ 

0  45+0-03 

^-0.03 

10.81? 

A  4+10 
-1.0 

46+10 
^  -10 

7  1+0-3 
'-4-0.2 

0  37+0-03 
'-0.03 

3-51! 

0  94+0-02 
44- ^^-0.02 

0-78!?:?! 

- 

0.88/48 

5 

g  04+0.05 
^•^^-0.06 

0  42+0-05 
-0.08 

10  7+0-4 
^^•'-0.6 

5-4!!:? 

52+^^ 

^  -17 

7  0+^-^ 

'  -O-0.5 

A  44+0.08 

^-^^-0.07 

3.6!?:? 

0  92+0-05 
44-^^-0.05 

0.94!?:?? 

- 

0.79/48 

6 

10.18t« 

-0.07 

11  3+0-5 

5.8!!:? 

51+^^ 

^4-13 

7  9+0.6 
'  -  -0.3 

0  31+0  02 
^-^4_o.o5 

3.0!?:? 

n  Q  9 +0.04 

44-^-+-0.04 

0-88!?:?? 

1  04+0-02 
4-44^-0.02 

1.09/75 

7 

g  40+0.03 
^-0.03 

0.36l« 

10.81? 

6  4+^’^ 

u.H-1  2 

61+1^ 

04-10 

7  7+0.7 
'-'-0.4 

0  30+0  04 
-0.04 

3  1+0-2 
'^•^-0.2 

0  94+0-02 
44- ^^-0.02 

0.84!?:?! 

1  05+^-^2 
•^^-0.02 

1.18/75 

8 

0-381“ 

10.21? 

5  9+^0 

61+10 

04-10 

7  7+0.4 
'-'-0.3 

n  99+0.04 
-0.04 

9  9+0.2 
•+-•+-0.2 

Q  93+0.02 

44-^-+-0.02 

0-78!?:?! 

- 

1.27/48 

9 

7  84+0-03 
-0.02 

0  50+0-03 
^•'^^-0.04 

11-51? 

5-1!!:? 

60+^0 

oO-io 

7  7+0.4 
'-'-0.3 

0  32+0-03 
-0.03 

9  9+0.2 
•+-•+-0.2 

Q  93+0.02 

44-^-+-0.02 

0  75+0-01 

44-  ''l-O.Ol 

0-96-X 
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Notes.  Fits  were  performed  with  the  CutoffPL  model  and  a  varying  number  of  CRSFs.  Uncertainties  and  upper  limits  are  at  the  90%  confidence 
level  for  one  interesting  parameter.  Epoch  for  data  grouping  (see  text  for  details).  Absorbed  luminosity  in  units  of  10^^  erg  s"^  Lx  covers  the 
energy  band  3-50  keV  and  was  calculated  using  a  distance  of  7  kpc  (Negueruela  &  Okazaki  2001).  In  units  of  photons  s"^  cm"^. 


can  be  used  to  describe  the  broad  band  spectra  of  accreting 
neutron  stars.  Due  to  the  early  successes  of  fitting  spectra  of 
4U  0115+634  with  the  NPEX  and  PLCUT  models,  these  continua 
have  been  extensively  used  when  modeling  the  outburst  behav¬ 
ior  of  this  pulsar.  Unfortunately,  however,  in  most  previous  anal¬ 
yses  -  including  some  of  our  own  -  no  attempt  was  made  at 
modeling  the  data  with  any  of  the  other  available  models. 

As  an  example,  modeling  the  2008  outburst  of  4U  01 15+634 
with  the  NPEX  model  (following  Nakajima  et  al.  2006,  fixing  T2 
to  2.0)  results  in  fits  that  are  of  only  slight  worse  quality  as  those 
using  the  CutoffPL  model  modified  by  the  10  keV  feature,  and 
we  can  recover  the  strong  variation  of  the  cyclotron  line  found 
in  earlier  analyses  (Fig.  6,  cyan  data  points)'^.  Depending  on 
the  choice  of  the  continuum,  we  therefore  find  a  fundamentally 


^  Modeling  the  spectra  with  the  PLCUT  model  results  in  an  unphysical 
break  around  the  energy  of  the  fundamental  line,  we  therefore  consider 
this  continuum  to  be  unsuitable  for  studying  the  cyclotron  line  behavior. 


different  behavior  of  the  fundamental  cyclotron  line,  especially 
the  centroid  energy. 

Given  that  the  fits  with  NPEX  and  CutoffPL  give  similarly 
good  results,  what  is  the  fundamental  difference  between  both 
modeling  approaches?  NPEX  models  generally  result  in  very 
broad  cyclotron  lines,  with  Wo  and  Wi  often  exceeding  5  keV 
(e.g.,  Nakajima  et  al.  2006,  runs  using  the  PLCUT  model  give 
similar  results,  see,  e.g.,  Li  et  al.  2012).  In  contrast,  using  the 
CutoffPL  continuum  with  the  10  keV  feature  results  in  widths 
of  typically  less  than  3  keV,  much  more  consistent  with  the  nar¬ 
row  shapes  of  the  residuals  shown  in  Fig.  2.  When  lines  are  as 
broad  as  in  the  NPEX  fits,  they  can  influence  the  continuum  fit. 
To  illustrate  how  strongly  these  features  actually  distort  the  con¬ 
tinuum,  Fig.  8  shows  the  continuum  shape  inferred  when  setting 
To,  Ti,  T2,  and  T3  to  zero.  The  residuals  for  the  CutoffPL  con¬ 
tinuum  with  the  10  keV  feature  (dark  blue)  show  sharp  absorp¬ 
tion  dips  at  the  cyclotron  line  energies  of  the  fundamental  and 
first  harmonic  line  energies.  For  the  NPEX  model  (green  line). 
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Fig.  5.  Left  temporal  evolution  and  right  luminosity  dependence  of  the 
continuum  parameters  (see  text  for  deflnitions)  and  10  keV  feature  for 
the  fits  with  the  CutoffPL  continuum.  The  color  gradient  from  reddish 
to  blueish  data  points  indicates  the  temporal  evolution  of  the  outburst. 
The  data  points  show  results  from  the  initial  fit  run,  where  a  variable 
number  of  cyclotron  lines  was  used.  Green  data  points  in  the  left  column 
correspond  to  the  second  fit  run  with  a  flxed  number  of  three  cyclotron 
lines  (only  in  the  left  column),  which  was  been  performed  as  a  consis¬ 
tency  check  (see  text  for  details).  Where  no  green  points  are  visible, 
both  fits  gave  the  same  results.  Lx  is  in  units  of  10^^  erg  s"^  Aq  and  Ape 
are  in  units  of  photons  s~^  cm“^,  while  Ape  is  in  multiples  of  10“^. 


on  the  other  hand,  it  is  obvious  that  rather  than  describing  nar¬ 
row  cyclotron  lines,  the  multiplicative  broad  line  model  strongly 
influences  the  continuum.  As  discussed  above,  our  spectral  fits 
show  that  during  the  brightest  phases  of  the  outburst  a  signifi¬ 
cant  hardening  of  the  continuum  emission  is  observed,  together 
with  a  change  of  the  exponential  cutoff.  While  this  is  also  seen 
in  the  variation  of  the  continuum  parameters,  it  is  very  likely  that 
the  luminosity  dependence  of  the  CRSF  in  the  NPEX  fits  is  due  to 
the  line  partially  modeling  this  behavior  of  the  continuum,  and 
not  a  real  physical  effect. 
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Fig.  6.  Luminosity  and  time  dependence  (inset)  of  the  energy  of  the 
fundamental  CRSF  using  the  CutoffPL  continuum  model.  The  color 
coding  is  the  same  as  in  Fig.  5.  Cyan  data  points  (diamonds)  show  the 
results  when  using  the  NPEX  continuum. 
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Using  the  results  from  the  CutoffPL  plus  10  keV  feature 
fits  for  the  line  energy  also  solves  another  problem  found  in  ear¬ 
lier  spectral  modeling.  Here,  using  NPEX  fits,  the  ratios  between 
the  fundamental  and  higher  cyclotron  lines  were  often  found  to 
deviate  significantly  from  integers.  While  slightly  non-integer 
ratios  are  expected  when  taking  relativistic  quantum  mechan¬ 
ics  into  account  (see  the  discussion  by  Pottschmidt  et  al.  2005), 
the  large  deviations  from  integer  multiples  seen  previously  (e.g., 
Tsygankov  et  al.  2007;  Nakajima  et  al.  2006;  Heindl  et  al.  1999; 
Santangelo  et  al.  1999)  require  rather  complex  model  assump¬ 
tions  which  introduce  second  order  effects  such  as  strong  vertical 
5-field  gradients,  crustal  field  structures  on  small  scales,  or  ther- 
momagnetic  effects  for  their  explanation  (see  Schonherr  et  al. 
2007,  and  references  therein  for  a  discussion).  Figure  9  shows 
that  we  can  recover  these  non-integer  ratios  in  our  NPEX  fits,  but 
when  using  the  CutoffPL  plus  10  keV  feature  model,  these  ra¬ 
tios  are  mostly  in  agreement  with  integer  values,  or  at  least  only 
slightly  higher,  as  expected  from  the  most  simple  models  for  cy¬ 
clotron  line  formation.  We  conclude  that  also  from  a  physical 
point  of  view,  the  CutoffPL  continuum  with  the  10  keV  feature 
gives  a  much  more  satisfactory  description  of  the  data. 

Finally,  we  check  the  dependence  of  the  line  position  on 
other  free  fit  parameters  when  using  the  CutoffPL  continuum 
modified  by  the  10  keV  feature.  Since  the  origin  of  the  10  keV 
feature  remains  still  an  open  questtion,  possible  relations  be¬ 
tween  the  CRSF  centroid  energy  and  the  parameters  of  the 
10  keV  feature  have  to  be  investigated,  and  other  continuum  pa¬ 
rameters  such  as  L’foid  or  F  could  influence  the  results  for  the  cy- 
lotron  line  parameters.  Figure  10  shows  the  behavior  of  Ax^  in 
the  vicinity  of  the  best  fit  values  for  three  epochs  representative 
for  the  early  phase,  peak,  and  late  phase  of  the  2008  outburst  for 
selected  combinations  of  important  fit  parameters.  These  contour 
plots  show  that  the  line  energy  is  not  affected  by  cross  correla¬ 
tions  with  other  parameters,  i.e.,  possible  variations  in  the  line 
energy  of  the  fundamental  cyclotron  line  are  not  due  to  variations 
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Fig.  7.  Time  (left)  and  luminosity  (right)  dependence  of  the  parameters 
of  the  CRSK  The  color  coding  is  the  same  as  in  Fig.  5. 


observed  in  the  parameters  of  the  10  keV  feature  or  the  contin¬ 
uum  parameters. 

In  contrast  to  the  stable  behavior  of  the  line  when  using 
CutoffPL  with  a  10  keV  feature,  Fig.  11  illustrates  the  depen¬ 
dence  of  the  CRSF  centroid  energy  on  the  continuum  parameters 
of  the  NPEX  model.  In  many  cases,  strong  correlations  between 
Eq  and  r  1  as  well  as  between  Eq  and  E’foid  are  present.  These  cor¬ 
relations  are  a  further  indicator  that  the  Lorentzian  component  is 
used  to  model  the  continuum  rather  than  a  cyclotron  line  when 
using  the  NPEX  model. 

5.  Summary  and  conclusion 

In  this  paper  we  have  presented  a  study  of  the  2008  outburst  of 
4U  0115+634  based  on  RXTE  and  INTEGRAL  data  (Eig.  1). 
We  reproduced  the  results  from  previous  work  (e.g.,  Nakajima 
et  al.  2006;  Tsygankov  et  al.  2007;  Li  et  al.  2012,  and  others), 
that  the  spectra  can  be  successfully  modeled  by,  e.g.,  the  NPEX 
model.  We  showed  that  for  these  fits  the  very  broad  absorption 
features,  which  are  thought  to  describe  the  CRSEs,  rather  model 
the  broadband  continuum.  Our  result  that  the  continuum  param¬ 
eters  are  strongly  variable  over  the  outburst  (Pig.  5),  could  be 
responsible  for  the  change  in  the  cyclotron  line  energy  in  these 
fits. 

We  have  shown  that  the  spectral  continuum  can  be  also  well 
described  with  a  model  introduced  by  Klochkov  et  al.  (2007)  for 
EXO  2030+375  and  Perrigno  et  al.  (2009)  for  4U  0115+634, 
i.e.,  the  CutoffPL,  modified  by  strong  Gaussian  emission  fea¬ 
ture  around  10  keV  and  several  (up  to  five)  cyclotron  lines 
(Table  2  and  Pig.  4).  Due  to  an  unphysical  break  close  to  the 
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Fig.  8.  a)  PCA,  HEXTE,  and  ISGRI  spectra  from  epoch  9  (red)  together 
with  the  best  fit  models  using  the  NPEX  (green)  and  CutoffPL  (blue) 
continua  after  setting  the  optical  depth  of  the  cyclotron  lines  to  zero,  b) 
Ratio  between  the  data  and  the  model. 


time  (MJD-54500) 

Fig.  9.  Ratios  of  the  cyclotron  line  energies  En  /Eq  with  respect  to  the 
centroid  energy  of  the  fundamental  line  against  time.  The  dark  green, 
dark  blue,  red  and  purple  data  points  (bars)  correspond  to  Ei  /Eq,  E2 /Eq, 
E^/Eq,  and  E^/Eq,  respectively,  using  the  CutoffPL  continuum  with  a 
10  keV  feature.  The  light  green  and  cyan  data  points  (filled  circles)  show 
the  results  for  Ei/Eq  and  E2/E0  using  NPEX.  The  dashed  lines  indicate 
integer  values  for  these  ratios. 


energy  of  the  fundamental  cyclotron  line,  the  PLCUT  model  does 
not  give  a  physical  description  of  the  continuum.  Using  the 
CutoffPL  together  with  the  10  keV  feature,  we  have  shown  that 
the  modeling  of  lines  reflects  the  theoretical  expectation  in  both 
the  line  shape  and  the  ratios  of  the  harmonics.  In  this  model,  the 
luminosity  dependency  of  the  centroid  energy  is  not  present.  As 
shown  also  by  Eerrigno  et  al.  (201 1),  the  line  centroid  energy  is 
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Fig.  10.  Confldence  contours  between  the  fundamental  CRSF  energy 
and  other  fit  parameters.  The  columns  show  the  results  for  epochs  2, 
6,  and  14,  respectively.  Contour  lines  correspond  to  the  68.3%,  90%, 
and  99%  confldence  level.  Color  indicates  with  respect  to  the  best 
fit  value,  with  the  color  scale  running  from  orange  (low  to  dark 
blue  (large 
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Fig.  11.  Confldence  contours  between  the  fundamental  cyclotron  line 
energy  and  continuum  parameters  when  using  the  NPEX  model.  The 
displayed  epochs,  colors  and  lines  are  the  same  as  in  Fig.  10. 


inferred  from  the  NPEX  and  PLCUT  models  was  almost  a  jump 
at  a  3-50  keV  luminosity  of  around  3x10^^  erg  s“^  and  strong 
hysteresis  effects  were  present  (Fig.  6).  Such  a  behavior  has  not 
been  seen  in  any  of  the  other  cyclotron  line  sources.  The  sec¬ 
ond  strongest  luminosity-dependent  cyclotron  line  variability  is 
observed  in  V0332-F53  (Mowlavi  et  al.  2006;  Tsygankov  et  al. 
2006;  Nakajima  et  al.  2010),  in  which  the  cyclotron  lines  are  so 
strong  that  the  change  can  be  seen  by  eye  even  in  the  raw  detec¬ 
tor  spectra.  While  the  boundaries  between  the  different  accretion 
regimes  estimated  by  Becker  et  al.  (2012)  might  therefore  have 
to  be  adjusted,  the  overall  principal  behavior  of  the  accretion 
column  discussed  by  these  authors  still  holds. 

With  the  recently  successful  launch  of  NuSTAR  (Harrison 
et  al.  2010),  a  low  background,  focusing  telescope  with  slightly 
higher  resolution  than  the  RXTE-HEXTE  or  the  Suzaku-HXD 
is  now  available.  With  its  higher  resolution  and  the  much  higher 
expected  signal  to  noise  ratio  in  the  spectra,  we  will  be  able, 
for  the  first  time,  to  resolve  cyclotron  lines  and  study  their 
luminosity-dependent  behavior  without  the  danger  of  a  system¬ 
atic  influence  of  the  choice  of  a  continuum. 


remarkably  stable  throughout  the  brightest  phase  of  the  outburst. 
Here  we  show  that  such  stability  is  present,  albeit  at  lower  sig¬ 
nificance,  also  during  the  early  and  late  phases.  In  this  model,  the 
continuum  variations  are  explained  by  a  slight  anticorrelation  of 
the  centroid  energy  of  the  Gaussian  feature  with  luminosity  and 
a  significant  variation  of  E. 

The  strong  systematic  influence  of  the  chosen  continuum 
on  the  cyclotron  line  behavior  illustrates  a  problem  that  is 
potentially  inherent  to  all  cyclotron  line  measurements.  We  em¬ 
phasize,  however,  that  this  does  not  mean  that  the  general  pic¬ 
ture  of  the  different  regimes  of  magnetized  neutron  star  accre¬ 
tion  and  general  luminosity  dependence  of  the  cyclotron  line 
as  outlined  by  Becker  et  al.  (2012,  and  references  therein)  is 
incorrect.  In  the  sample  of  sources  discussed  by  Becker  et  al. 
(2012),  4U  01 15-f634  was  an  outlier.  The  change  of  line  energy 
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